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Thermoelectric transport in strongly spin-orbit coupled two-dimensional Rashba system is studied 
using the exact solution of the linearized Boltzmann equation. Some unusual transport behaviors 
are revealed. We show that the electrical conductivity takes a Drude form when the Fermi energy 
Ef is above the band crossing point, but a non-Drude form which is a quadratic function of Ef when 
Ep lies below the band crossing point. The Mott relation breaks down when Ep lies in the vicinity 
of the band crossing point. It is shown that the thermopower and thermoelectric figure of merit are 
strongly enhanced when Ep downs below the band crossing point. This enhancement is attributed 
to not only the one-dimensional-like density of state but also the unconventional intraband elastic 
scattering below the band crossing point. The differences between these results and those obtained 
by the relaxation time approximation are discussed in detail. 
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I. INTRODUCTION 

In two-dimensional electron system (2DES) with 
Rashba spin-orbit coupling (SOC), when the energy 
downs below the band crossing point (BCP), a band val¬ 
ley emerges (Fig. 1) with different topology of Fermi 
surfaces (FS) from that when the energy is above the 
BCP— 3 . In the band valley the dispersion curve is 
not monotonic and the density of state (DOS) is one¬ 
dimensional (lD)-like. In the case of strong SOC, the 
Fermi level can lie in the vicinity of or even below the 
BCP, and the valley structure can survive the weak dis¬ 
order and thermal smearing at low temperatures. There¬ 
fore, in this case the nontrivial topology of FS in the band 
valley may affect the transport properties significantly. 

The strongly spin-obit coupled 2DES, formed at the 
Te-terminated surface in layered polar semiconductors 
BiTeX (X=C1, Br, I), has been discovered recently 
by ARPES measurements in agreement with ab-initio 
calculations^—. In such 2DES the giant Rashba SOC 
coefficient is one order of magnitude higher than that in 
conventional III-V semiconductor heterostructures*. In 
addition, the BiTeX quantum well^ may be another can¬ 
didate to realize strongly spin-orbit coupled 2DES. Very 
recently the first-principle band structure calculation has 
suggested the strain engineering of heavy-metal film on 
layered large-gap semiconductor substratold. as a promis¬ 
ing way to form 2DES with strong Rashba SOC, e.g., Au 
single layer on strained InSe(OOOl). In aforementioned 
2DES, the large Rashba spin-splitting provides a chance 
to study the unconventional transport properties induced 
by the band valley and band crossing. 

The effects of the band valley and band crossing on spin 
transport and superconducting electronics have received 
much theoretical attention, e.g., the non-Dyakonov-Perel 
spin relaxation behavior—, the non-Edelstein electric- 
field induced spin polarization— the enhanced spin- 
orbit torque efficiency^, the enhanced superconducting 


instabilities^ and the specular Andreev reflection in the 
interface of a superconductor and a 2DES with strong 
Rashba SOC— In addition, spin-related thermoelectric 
conversion in systems with strong Rashba SOC is gath¬ 
ering increasing attention, which is not only essential for 
exploring spintronics devices^ but also important for de¬ 
velopments of spin caloritronics— For Rashba 2DES, 
based on the relaxation time approximation (RTA)£ i 10 i 18 
in the semiclassical Boltzmann equation (SBE) approach, 
it has been suggested recently that the dimensional re¬ 
duction of the electronic structure from 2D to ID can 
result in enhancements of the diffusive thermopower and 
thermoelectric figure of merits ZT = [a/a) 2 aT / n. 
Here er, a , k and T denote the electrical conductivity, 
Peltier coefficient, thermal conductivity, absolute tem¬ 
perature, respectively. In the RTA, the enhancement of 
thermopower was attributed solely to the ID-like DOS 
below the BCP—. 

However, so far, no fully satisfactory theoretical study 
on thermoelectric transport exists for the case that Ep 
lies in the vicinity of or below the band crossing point. 
This is because that the RTA is inappropriate for Rashba 
2DES at low temperatures when the electron-impurity 
scattering dominates, in the case of strong SOC—. For 
Fermi energies above the BCP, the RTA can not han¬ 
dle the difference in the relative importance between the 
interband and intraband elastic scatterings. This differ¬ 
ence is significant when Ep lies near the BCP, so the 
RTA is unsuitable for this case. While for Fermi ener¬ 
gies below the BCP, the nontrivial FS topology induces 
not only the ID-like DOS, but also nonconventional in¬ 
traband scatterin@I4(inter-branch and intra-branch scat¬ 
terings, the two branches are denoted in Fig. 1) which is 
also beyond the scope of RTA. It has been shown thal—, 
when Ep is below the BCP, the nonequilibrium spin po¬ 
larization calculated by an exact transport time solution 
of the SBE is quite different from the result obtained by 
the constant RTA—. This motivates us to employ the 
exact solution— of the SBE to systematically study the 
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effects of strong Rashba SOC on the spin-related ther¬ 
moelectric transport in 2DES at low temperatures. 

In this paper the exact solution of the SBE is employed 
to calculate thermoelectric transport coefficients and the 
figure of merit in Rashba 2DES. This solution is suitable 
when electron-impurity scatterings dominate. We show 
that the electrical conductivity takes a Drude form when 
Ep is above the BCP, but a non-Drude form which is a 
quadratic function of Ep when Ep is below the BCP. We 
found that, the .Ep-dependence of the Peltier coefficient 
is not monotonic and the Mott relation^ breaks down in 
the vicinity of the BCP. It is shown that the thermopower 
and thermoelectric figure of merit have strong enhance¬ 
ments when Ef is tuned below the BCP. This enhance¬ 
ment in the thermoelectric performance is a combined 
result of the ID-like DOS and the unconventional intra¬ 
band scattering induced by the nontrivial FS topology in 
the band valley. 

The paper is organized as follows. We present the char¬ 
acteristic properties of Rashba 2DES in Sec. II. The ther¬ 
moelectric transport coefficients are calculated in Sec. 
III. In Sec. IV, the thermoelectric figure of merit is pre¬ 
sented. The conclusions of the present paper are given 
in Sec. V. 


II. CHARACTERISTIC PROPERTIES OF 
RASHBA SPINTRONIC SYSTEMS 

A. The DOS and the topological features of FS 

We consider a Rashba 2DES with spin-independent 
disorder 

iJ= ^ + f^' (pX " ) + F(r) ’ (1) 

where p = fik is the momentum of the electron, m is 
the effective mass, a = (a x ,a y ,a z ) are the Pauli ma¬ 
trices, /3 the Rashba coefficient. The disorder potential 
V(r) = EiWr-Ri) is produced by randomly dis¬ 
tributed identical 5-scatters at and is modeled by 
the standard Gaussian disorder average (|Vk'k| 2 ) = 

n-imYo where rq m is the impurity concentration, Vq is the 
strength of the disorder potential, Vk' k is the orbital dis¬ 
order matrix element and {--) dis the disorder average. The 
inner eigenstates and eigenenergies of the clean system 
read |zt Ak ) = ^ (1, -*A exp (i</>)) T and E xk = ^ +X/3k , 
respectively. Here A = ± and <j> is the polar angle of k. 

The DOS at energy E > 0 is given by 7V> ( E ) = 
£ a 7V a (E) = 2No with N\ (E) = N 0 k ^l kR the 
DOS in the A band. Here kp = m jp, N 0 = 
k\ ( E ) = —A kp + -g E z r + 2 EpE is the wave number 
corresponding to a given energy E > 0 in the A band, 

Ep = m is the ”Rashba energy”. 

Below the BCP there is a valley structure in the lower 
Rashba band, with the bottom located at kp and the 
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FIG. 1. Band structure of the Rashba 2DES. (a) Dispersion 
curve. The energy of the bottom of the dispersion curve is 
E- (Icr) = —\Er with kp = m-P?. At a given energy E > 0, 
the wave number in ± band is defined as k± (E). For — ^Ep < 
E < 0, there are two monotonic regions on E — k curve: the 
one from k = 0 to kp is marked by the branch —2, whereas the 
other from k = kp to 2 kp marked by branch —1. The wave 
number ( E ) represents the wave number in the — v branch 
at given E, where v = 1,2. (b) Constant-energy circles for 
E > 0. (c) Constant-energy circles for —^Ep < E < 0. The 
arrows in (b) and (c) represent the directions of the group 
velocity. 

minimal energy E_ (kp ) = —^Ep. At a given energy 
E in the band valley above the bottom, there are two 
wave numbers fc_„ (E) = kp + (—l)" -1 ^ \JE^ + 2 EpE 
where v = 1,2 denote the two monotonic branches as 
shown in Fig. 1. It is worth noticing that the group 
velocity in the —2 branch is anti-parallel to the mo¬ 
mentum. The DOS in this regime has a ID behav¬ 
ior N k (E) = E 2 =i N-AE) = ™ 0 Ve * + « ErE with 

N-AE) = N 0 lk _ k - { ^% R \ . 

B. The response to the external electric field and 
temperature gradient 

In the linear response regime, within the semiclassical 
Boltzmann framework the electric current density j e and 
heat current density j h are given by 

j e = e E/‘ v * = aE *+a(-VT) (2) 

l 

and 

j h = Z!/i v «(E|-/i) = Ta E *+«(- v r). (3) 

i 

Here l is the eigenstate index, v; is the group velocity of 
state l, fi denotes the semiclassical distribution function 
(DF) for the electron wave packets, E* = E — AV/r the 
gradient of the electrochemical potential, VT the temper¬ 
ature gradient, fi = f° (E;) + gi where /° is the equilib¬ 
rium Fermi-Dirac DF and gi the out-of-equilibrium devi¬ 
ation linear in the generalized driven force. The system is 













3 


time-reversal invariant, thus the Hall transport is absent 
and cr, a, k are all numbers. To calculate these thermo¬ 
electric transport coefficients in low temperature cases 
where the static impurity scattering dominates, we em¬ 
ploy the SBE^. In the presence of weak uniform electric 
field and gradients of chemical potential and tempera¬ 
ture, the linearized SBE suitable for the present system 
in nonequilibrium stationary state takes the formlS^: 

d 

Fl '' ri 7rE = -'52 UJl, ’ l { 9l ~ 91 '}' ( 4 ) 

^ 7 / 


Here F( = eE* + E ‘~ M (—VT) is the generalized force 
acting on state l. uji^i is the elastic scattering rate 
from eigenstate l to l', which can be obtained by 
the golden rule in quantum mechanical scattering the¬ 
ory. For the present system, the Born approxima¬ 
tion in the lowest order is sufficient^ 21 *, i.e., = 

f (|WV \{uv\ui)\ 2 5 (El - E u ). 

The SBE in isotropic 2D Rashba system can be solved 
conveniently using energy E and polar angle 4> and band 
index A as variables, i.e., I = (E, A, 4 1 ). While the val¬ 
ley region of the lower band is worth noticing due to the 
non-monotonic band dispersion. In the valley region, the 
branch index —v introduced above is needed to denote 
the eigenstate, i.e., I = (E,—v,(j)). Then the exact solu¬ 
tion of Eq. 0 can be obtained, which we refer to our 
previous workM: 

9x (E) = (~d E f°) F e ■ k ^lr (5) 

v 7 m 

when E > 0, and 


9-u (E) = (-d E f°) F E Hk ~ AE) T (-1) 

m 


E\ + 2 E r E 

El 


( 6 ) 

when E_ (k R ) < E < 0. Here F^ represents the 
generalized force acting on electrons with energy E, 

r = ^ 2;rll ‘” 1 o N ° ^ i s the ordinary momentum relax¬ 


ation time. At E = 0, g + (E —> 0 + ) = g- 2 (E — > 0 ) = 
0, 5 _i(E-s>0-) = g-{E^ 0+). ' The DOS at 

(E = 0, k = 0) vanishes, so this point does not contribute 
to transport quantities in Eq. 0 and 0, and only 
the outer constant-energy circle (E = 0, k = 2 k R ) con¬ 
tributes at E = 0. 

We note that Eq. 0 takes into account the intraband 
(A —> A) and interband (A —> —A) elastic scatterings. 
While Eq. 0 takes into account the unconventional 
intraband scattering in the band valley, i.e., the intra¬ 
branch (—v —> —v) and inter-branch (—u —> —(3 — v)) 
elastic scatterings. 


C. The chemical potential at low temperatures 

The relation between E e and g is needed for investi¬ 
gating the Ep -dependence of Peltier coefficient and ther¬ 
mal conductivity. It can be obtained by considering the 


electron density as follows. At finite temperatures, the 
electron density can be calculated by n e = nf + nf, 
where we define n> and nf as 


= 


pOO 

/ dEf° (E) 7V> (E ), 

Jo 

L 


(7) 


E-(k R ) 


dEf° (E) iV< (E). 


By defining two functions: 

l>E 

Q> (E) = / dE’ 7V> (E') = 2N 0 E (8) 

Jo 

for E > 0 and 


Q< (E) = j dE'N < (E 1 ) = 2N 0 E r f 1 


1 + 2 


E 

~E~r 


, (9) 

for E < 0, nf and nf can be written in the following 
forms via integration by parts: 


n> =dE(~ypjQ > (E), ( 10 ) 

n< = Q< (E_ (k R )) - £ ^ dE £j£) Q< ( E ). 


Setting 


E - g _ g _ E F 
k B T ~ X: k B T~ lj k B T ~ 


( 11 ) 


in Eq. m and restricting to not too low chemical po¬ 
tential ^ 1, the total electron density is found 

as 


n e 

2NqE r 



ti) 


k B T 

Er 


+ 1 



l + 2(:c 


t\) 


k B T 
Er ' 


( 12 ) 


Here the condition M f \ can kg realized at low 

temperatures due to the giant Rashba SOC, e.g., in BiTel 
surface stated (k R ) ~ —90meF, if g = \E_ (k R ) we 
have ~ 450AT. At zero-temperature g = Ep , 

when Ep > 0 it has been obtained that^ 2 n e (Ep > 0) = 
2-ZVq (Ep + Er). While when Ep < 0 we obtain 


n e (E f < 0) = f 

Jk 


k-i(E F ) 


k-2(EF) 


kdk 

27r 


= 2 Nr 


El 


2 ErE f . 


(13) 

In this case the Fermi level intersects only the lower band 
and only the annulus lying between the two Fermi circles 
of radii /c_i ,_2 (Ef) is filled. This nontrivial topology of 
Fermi surfaces in the band valley has been highlighted in 
previous researches^* 3 . Substituting the electron density 
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E f ( 1CT 2 E r ) 


FIG. 2. The difference between the chemical potential and the 
Fermi energy. The temperature is fixed to 0.02 Er/Icb- The 
numerical results describe exactly the behavior near Ef = 0, 
which is beyond the scope of the Sommerfeld expansion. 


into Eq. m, we obtain the Fermi energy dependence of 
chemical potential at a given temperature, as presented 
in Fig. 2. 


When Ef > 0.1 Ep, /i — Ep = 0, same as the analytic 
result based on the Sommerfeld expansion^ in the case 
of Ef ksT. When Ep < — 0.08Er, numerical results 
can be well fitted by the formula /r — Ep = if e r + 2 E f 
obtained by the Sommerfeld expansion in the case of 
—Ep kpT and Ep + \Ep kpT (this latter con¬ 
dition ensures that the band valley structure and Fermi 
surface are not smeared out by thermal broadening and 
can be satisfied at low temperatures with giant Rashba 
SOC). It is well-known that at low temperatures in three- 
dimensional (3D) parabolic 2DES p, < Ep, in 2D /i = Ep 
and in ID p, > Ep , due to different energy dependencies 
of the DOSi^. Here the above two limiting cases corre¬ 
spond to 2D and ID cases, due to the 2D- and ID-like 
DOS for energies above and below the BCP, respectively. 


When Ep is in the intermediate region — 0.08 Er < 
Ef < 0.1 Ep, the Sommerfeld expansion is not valid and 
numerical results clearly show continuous transition be¬ 
tween above two limiting cases. This transition from 2D 
to ID is resulted by the Fermi surfaces topology change 
in Rashba model. 


III. THERMOELECTRIC RESPONSE 
COEFFICIENTS 


A. Drude-like and non-Drude forms of electrical 
conductivities 


In what follows we assume that the generalized force 
is applied in x direction. When both bands are par¬ 
tially occupied, the zero-temperature electrical conduc¬ 
tivity is calculated by substituting the group velocity 
v (2£, A, (f>) = jvf(g) an d th e DF Eq. 10 into Eq. 


©■ The result is 


a (Ep > 0) 


e 2 27T (Ep + Er) t 
2'k 2 h h 


(14) 


where and below we use the notation a ( Ep) to 
represent the zero-temperature electrical conductivity 
a (T = 0, Ep) for brevity. 

This result still has the usual form of Drude formula 
(7 = n e e 2 r/m. It has been obtained in some earlier works 
based on the Green’s function calculation in the ladder 
approximation^ or the same exact solution to the SBE^ 
as our Eq. 0 when Ep > 0. In this case it looks as if 
the electric current were generated by charge carriers of 
one type with density n e and mobility er /m. 

For Fermi energies below the BCP, substitut¬ 
ing g_v and the group velocity v (E, —u, </>) = 

(—l)" -1 N N ° E j i n t° Eq. 0, we obtain the zero- 

temperature electrical conductivity 


m ^ rw e 2 27 t (Ep + E R ) t E^ + 2 E R E F 

a{EF - Q)= 2^h - h - & R -• 

(15) 

It has a quadratic dependence on the Fermi energy and 
does not take the form of Drude conductivity, different 
from Eq. ©. Since the Fermi surfaces topology in the 
band valley differs from that above the BCP, the behav¬ 
iors of electrical conductivities are different between the 
two regions. 


B. The Peltier coefficient 


The Peltier coefficient and the electrical conductivity 
are connected by (details in the Appendix) 



dE 


’E_(k R ) 



E — g, 
T 


a (E) ! 


( 16 ) 


where cr (E) represents the zero-temperature electrical 
conductivity with Fermi energy E. Usually, this energy 
integration is worked out by performing the Sommerfeld 
expansion^. When the Sommerfeld expansion is valid, 
the Mott relation holds. The former demands that cr (E) 
is continuously differentiable at the chemical potential. 
However, according to Eqs. m and (ED, cr (E) takes 
different forms between the two sides of the BCP. There¬ 
fore 1 F is not continuous and the Sommerfeld ex¬ 

pansion is not valid at /i = 0, thus the Mott relation fails 
for chemical potentials near the BCP. 

Substituting Eqs. (flTl) and (fl5l) into Eq. (fTCl) . we 
obtain the Peltier coefficient as 


= 7r 2 fc|T cr(0) r _ &(*!) -ha (U) 

3e Ep, \ 7T 2 

^~2 c(U) + 2*1 - 6(ii)^ +t\a (h) 


k B T 

Ep 


( 17 ) 
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FIG. 3. (a) The Fermi energy dependence of the Peltier Coefficient, (b) The comparison between the numerical results and 
the Mott relation. The deviation from the Mott relation is significant when the Fermi energy lies near the BCP Ef = 0. In 
(a) and (b), the temperature is fixed to 0.02E B /kB . The relation between the chemical potential and Fermi energy shown in 
figure 2 has been taken into account 


where we define 

a(ti) 

b(ti) 

c{t l) 



(18) 


Eq. m can describe the behavior of a for Fermi energies 
near the BCP. Because the Sommerfeld expansion is un¬ 
suitable in the band crossing region, we have to perform 
numerical calculations. 

For temperature at 0.02 Ep/k B , the Fir-dependence of 
Peltier coefficient is shown in Fig. 3(a). We find that 
when Ep/E r > 0.1, a is almost constant with respect to 
Ep , which is consistent with the Mott relation (Eq. (1191) 1: 
when Ep/E r < —0.08, a has a nearly linear dependence 
on Ep , which is also consistent with the Mott relation 
(Eq. (1201) 1. Between above two regions, a non-monotonic 
Fermi energy dependence of a in the band crossing region 
is found. 

We use the notation a Mott to denote the Peltier 
coefficient obtained by the Mott relation a Mott = 

n 2 k 2 B Tda(E), 

3e dE \ E—Ef • 


a Mott 


n 2 kgT a (0) 
3e Er 


E f > 0, 


(19) 


and 


a Mott 


K 2 k% T 
3e 



1 + 


4 Ep 
3 F^ 


, Ep < 0. 


( 20 ) 


Combining Eq. (1T71) with Eqs. (fTPl) and (1201) . we plot 
a/a Mott in Fig. 3(b). When -0.08 Er < E f < 0.1 Er, 
our results show deviations from the Mott relation. 

The Mott relation fails near the BCP due to the fact 
that the electrical conductivity takes different forms on 


the two sides of the BCP. Therefore the deviation from 
Mott relation can be regarded as a consequence of the 
topological change of FS varying across the BCP. 


C. The thermal conductivity 


The thermal current response to the temperature gra¬ 
dient can be obtained by (details in the Appendix) 




dE - 


E-(k R ) 


a/o 

dE 


E — fi 
k B T 


'(E)- 
( 21 ) 

Substituting Eqs. 1121) and (fl5l) into Eq. (12T1) . we obtain 
2c(ti) -Mi (tt 2 - 2 b(h)) k B T 


= 1 - 


Lqct (0) T ~ tt 2 /3 Er 

~d(t i)+ 2f lC (fi)-M? (x - 6 (*i)) / k B T\ 

^/6 


lul 

15 


V Er J ’ 
( 22 ) 


where Lq = 'f L ) 2 is the free-electron Lorentz num¬ 
ber, and d(ti) = x 4 . This is our main 

result for k. 

The Fp-dependence of k is given in Fig. 4(a). It shows 
monotonic dependence on the Fermi energy of k. 

We compare our numerical results with that given by 
the Wiedemann-Franz law in Fig. 4(b). Here k wf de¬ 
notes the thermal conductivity based on the Wiedemann- 
Franz law k wf = LTa(fj) ~ LTa (Ep). The deviation 


from Wiedemann-Franz law is very slight: 


< 


6% in the whole regime we investigate here. Hence 
the Wiedemann-Franz law holds quite well. The differ¬ 
ence between the validity of the Mott relation and the 
Wiedemann-Fran law in the band crossing region can 
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FIG. 4. (a) The Fermi energy dependence of the thermal conductivity, (b) The comparison between the numerical results and 
the Wiedemann-Franz law. The deviation is very slight, about < 6% near the BCP. In (a) and (b), the temperature is fixed to 
0.02En/kB- 


be understood as follows. The Mott relation connects 
the Peltier coefficient with the energy derivative of the 
electrical conductivity, while the Wiedemann-Franz law 
connects the thermal conductivity with the electrical con¬ 
ductivity itself. At the BCP, the electrical conductivity is 
continuous, but its energy derivative is not. The Peltier 
coefficient based on the Mott relation is not continuous 
across BCP. Therefore, in the vicinity of the BCP the 
Mott relation breaks down, while the Wiedemann-Franz 
law is satisfied quite well. 


D. Comparison with RTA results 


Eqs. (fl4l) and (fl5l) are different from the electrical 
conductivities obtained by employing the constant RTA-: 


a RTA (Ep > 0) = 


e 2 2 tt (E f + \E r ) t 
2n 2 h h 


a RTA (Ep < 0) 


e 2 r 
27 rh 2 


E\ 


2 E F E F . 


(23) 


The RTA result neither takes the Drude form for E F > 0 
nor has a polynomial dependence on Ep /Er for E F < 0. 
When E F > 0, a RTA tends to Eq. ([Till for weak SOC 
Ep Eft. For Fermi energies near the BCP, the elec¬ 
trical conductivity obtained by the constant RTA signif¬ 
icantly differs from Eqs. (fTil) and (fl5l) : cr (Ep = 0) = 
n e (Ep = 0)e 2 T/m = 2a RTA (Ep=0). For Ep be¬ 
low the BCP, as long as the band valley structure can 
survive the thermal smearing and disorder broadening: 
Er kpT, H/t, the difference between Eq. (fl5l) and 

<J RTA can n0 ^ k e ig norec J 

These differences between our results and the RTA re¬ 
sults can be understood as follows. When Ep lies high 
above the BCP Ep Er, N + (Ep) ~ N- (Ep), the 
intraband and interband scattering events are of equal 
importance for both the inner (+) and outer (-) Fermi cir¬ 
cles. Only in this case the RTA works well. When Ep lies 
near the BCP, N + (Ep) tends to zero. Thus, for Fermi 


electrons on the inner Fermi circle the interband scat¬ 
tering events dominate over the intraband scattering^; 
while for Fermi electrons on the outer Fermi circle the 
intraband scattering events dominate over the interband 
scattering. The difference in the relative importance be¬ 
tween the intraband and interband scattering events as 
well as its change with varying Ep can not be described 
by the RTA. When Ep lies below the BCP, the non- 
Drude form of <r RTA (Ep < 0) is caused only by the 1D- 
like DOS below the BCP, while that of Eq. (fl5l) based on 
the exact solution of SBE relies on not only the ID-like 
DOS but also the unconventional intraband scattering 
(inter-branch and intra-branch scatterings) induced by 
the nontrivial FS topology. 

Now we examine the thermopower in the ’’Mott rela¬ 
tion regimes”: S = a Mott /a. When Ep kpT 

s _ ^klT 1 _ p 2 k 2 B T 2N 0 
3e Ep + Eft 3e n e 

and when — Ep kpT 


S = 


n 2 k 2 B T 3 
3e Eft 


i _i_ 4 Ef 
E r 



(25) 


The thermopower is enhanced in the band valley, similar 
to the thermoelectric figure of merit which will be dis¬ 
cussed in the next section (Eqs. (071) 1. Eqs. (l24l) and 
(ESI are different from those obtained by the constant 
RTAi£: 


S RTA (Ep » k B T) 
S RTA (- Ep > k B T) 


n 2 k 2 B T 1 
3e E f + i Ep ’ 
17 r 2 k B T 1 

2 3e Ep + \Er 


(26) 


Based on the RTA result Eq. (l26l) . it had been concluded 
that below the BCP the enhancement of thermopower in 
a Rashba 2DES compared to a parabolic 2DES with the 
same electron density is caused solely by the much lower 
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Ep. And this lower Fermi energy is a consequence of the 
ID-like DOS^. However, Eq. (05l) obtained using the ex¬ 
act solution of the SBE shows that the enhancement of 
thermopower is a combined result of the ID-like DOS and 
the unconventional inter-branch and intra-branch scat¬ 
terings in the band valley. 

In conventional semiconductor asymmetric quantum- 
wells the difference between Eq. d and the 1st equa¬ 
tion of Eq. (l23l) is negligible due to the weak Rashba 
SOC Ep Ep, so does the difference between Eq. 
and the 1st equation of Eq. (l26l) . In Rashba semi¬ 
conductors BiTeX (X=C1, Br, I), despite that the re¬ 
ported Fermi levels in existing experiments are still in 
3D bulk conduction band (BCB), the validity of con¬ 
stant RTA analysis of electrical conductivity has been 
questioned when Ep lies near the BCP of BCB^. With 
further studies on systematic doping in these Rashba 
semiconductor s 24 ! 25 , it is promising that the Fermi level 
can be tuned into the bulk band gap. In addition, the 
BiTeX quantum welllS is another possible candidate to 
realize strongly spin-orbit coupled 2DES. Very recently 
the first-principle calculation has suggested the forma¬ 
tion of 2DES with large Rashba SOC by strain engi¬ 
neered growth of a Au single layer on the layered large 
band-gap semiconductor InSe(OOOl) substrate^. Future 
studies of the strain engineering of heavy-metal film on 
layered large-gap semiconductor substrate may also re¬ 
alize 2DES with stronger Rashba SOC. In these systems 
the transport properties of the 2DES with strong Rashba 
SOC can be detected experimentally and our theoretical 
results can be tested. For the experimental measure¬ 
ments of the low-temperature diffusive thermopower, the 
hot-electron thermocouple technique could be applied, 
which is much less sensitive to phonon-drag effects than 
conventional methods^. 


IV. THE THERMOELECTRIC FIGURE OF 
MERIT AND THE ENHANCEMENT BELOW 
THE BCP 


The performance of a thermoelectric material is deter¬ 
mined by the figure of merit ZT = (a/a) 2 aT/ k. 

For the case that the Fermi energy lies near the BCP, 
the ZT is shown in Fig. 5. It shows that when the Fermi 
energy is tuned across the BCP from Ep = 0.1 Ep down 
to Ep = —0.08 Ep, the figure of merit acquires a strong 
enhancement. 

When Ep/Ep > 0.1 or Ep/Ep < —0.08, the Mott 
relation holds very well, thus ZT can be calculated using 
a Mott , yielding ZT = (a Mott /a) 2 / L, i.e., 


ZT 


3 



Ep > O.IEr 


, Ep < —0.08 Ep 


( 27 ) 



FIG. 5. The Fermi energy dependence of the figure of merit 
ZT when the temperature is fixed to 0.02 Er/Icb- 


In both energy intervals, ZT monotonically increases 
with decreasing Fermi energy, similar to the case when 
—0.08 < Ep/Ep < 0.1 shown in Fig. 5. 

This enhancement in ZT when Ep is tuned below the 
BCP is directly related to the enhanced Peltier coefficient 
and the decreased electrical conductivity when Ep varies 
below the BCP, as shown in Fig. 3(a) and Eqs. (fldl) and 
(Ha. Because the Peltier coefficient is connected with 
the electrical conductivity , the enhancement in ZT can 
be attributed to the nontrivial topology of Fermi surfaces 
in the band valley regime. 


V. CONCLUSIONS 

We have calculated thermoelectric transport coeffi¬ 
cients and the figure of merit in strongly spin-orbit cou¬ 
pled Rashba 2DES with spin independent disorder using 
the exact solution of the linearized Boltzmann equation. 
At low temperatures the static impurity scattering dom¬ 
inates, it is shown that the electrical conductivity takes 
a Drude form when the Fermi energy Ep is above the 
band crossing point, but a non-Drude form which is a 
quadratic function of Ep for Fermi energies below the 
band crossing point. This is attributed to the different 
topologies of Fermi surfaces on the two sides of the band 
crossing point. For Fermi energies near the band crossing 
point, the Ep’-dependence of the Peltier coefficient is not 
monotonic, and the Mott relation breaks down. While 
the thermal conductivity is monotonically increasing as a 
function of Ep and the Wiedemann-Franz law holds quite 
well. The thermopower and figure of merit are strongly 
enhanced when Ep downs below the band crossing point. 
This enhancement is caused not only by the ID-like den¬ 
sity of state but also by the unconventional intraband 
elastic scattering below the band crossing point. 

Our results differ from previous ones obtained by the 
relaxation time approximation, especially for Fermi ener¬ 
gies in the vicinity of and below the band crossing point 
in systems with strong Rashba spin splitting. For Fermi 








energies above the band crossing point, our results can 
handel the difference in the relative importance between 
the interband and intraband elastic scattering events, in 
contrast to the relaxation time approximation. This dif¬ 
ference is significant when Ep lies in the vicinity of the 
band crossing point. For Fermi energies below the band 
crossing point, our results take into account the uncon¬ 
ventional intraband scattering induced by the nontrivial 
FS topology, which can not be described by the relax¬ 
ation time approximation. 

Our theoretical results may be tested in strongly spin- 
orbit coupled 2DES, e.g., the surface state of polar semi¬ 
conductors BiTeX (X=C1, Br, I) and BiTeX quantum 
wells, as well as the 2DES formed by the strain engineer¬ 
ing of heavy-metal film on layered large-gap semiconduc¬ 
tor substrate^. 
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Appendix A: Derivations of the analytic formulas 
for the Peltier coefficient and thermal conductivity 

In order to establish the connection between the elec¬ 
trical conductivity and Peltier coefficient and thermal 
conductivity, only the original transport time form of gi 
is needed, i.e., 


The thermal conductivity is found as k = k > + k < with 


df \ (E — n) 
dE 


=/Sfn- 

xJ2Nx(E)v 2 x (E,\,0)t x (E) (A7) 
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gx(E)= (-^pjF E .yr(E,X,ct>)Tx(E), (Al) 

and 

9-v (E) = F e ' v(E,-v,(/>)t- v {E ), (A2) 


Therefore the Peltier coefficient can be expressed as 


and 




E — g 
T 


a> (E) 


E — g 
T 


a(E) 


(A9) 


(In this model t x (E) = r N ^' > an d T -i/ (E) = 

r N ~^ E ^ Er + pf R E ■ which has been included in and can 
be read out from the specific form of nonequilibrium DF 
in the paper. However, for the purpose in this sup¬ 
plementary material, this specific form is not needed.) 
Substituting Eqs. JAD and (IA2I) into the two linear re¬ 
sponse equations, since the generalized force is applied 
in x direction, the the electrical conductivity is given by 
a = a > + a < , where 



xJ2Nx(E)v 2 x (E,X,cI>)tx(E) (A3) 
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so that 
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Similarly, the themal conductivity is connected with the 
zero-temperature electrical conductivity as 
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